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ABSTRACT 

We report the spectral analysis of the X-ray pulsar LMC X-4 in its high state 
out of eclipse observed by BeppoSAX. During this observation no coherent pulsations 
are detected. The primary continuum is well described by a power law with a high 
energy cutoff (EcutoS ~ -E'foid ~ 18 keV). The addition of a cyclotron absorption line at 
~ 100 keV improves the fit significantly. The inferred magnetic moment is 1.1 x 10^^ 
Gauss cm^, in agreement with the value estimated assuming that the neutron star is 
at the spin equilibrium, as it has been proposed for this source. The remaining excess 
at low energies can be fitted by a Comptonization of soft photons by moderately hot 
electrons (kT ~ 0.9 keV), with an optical depth r ~ 16. The seed photons for this 
Comptonization are consistent with black body emission from the accretion disk at the 
magnetospheric radius. Another possibility is to fit the soft excess with black body and 
thermal bremsstrahlung. In this case the black body would originate from cold plasma 
at the magnetosphere while the bremsstrahlung component may be produced by the 
strong stellar wind from the companion star, ionized by the X-ray emission from the 
pulsar. 

Subject headings: stars: individual: LMC X-4 — stars: magnetic fields — stars: neu- 
tron — X-rays: stars 

1. Introduction 

LMC X-4 is an eclipsing high-mass X-ray binary pulsar in the Large Magellanic Cloud. Its 
optical counterpart was identified with a 20 Mq 07 III-V star (Sandaleak &; Philip 1977). The 
orbital period of the system is 1.4 days (Li, Rappaport & Epstein 1978; White 1978). The X-ray 
intensity varies by a factor ~ 60 between high and low states with a cycle of 30.3 days (Lang et al. 
1981). This long-term variation is attributed to a periodic blockage of the direct X-ray beam by a 
precessing accretion disk, which is tilted with respect to the orbital plane of the binary (Lang et 
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al. 1981; Ilovaisky et al 1984, Priedhorsky & Holt 1987; Woo, Clark & Levine 1995, and references 
therein). Flaring episodes occur about once per day, lasting from ~ 20 s to 45 minutes, during 
which the intensity increases by factors up to ~ 20 (Epstein et al. 1977; White 1978; Skinner et al. 
1980; Kelley et al. 1983; Pietsch et al 1985; Dennerl 1989; Levine et al. 1991). 

Coherent pulsations with a period of 13.5 s were discovered by Kelley et al. (1983) during flares 
events. Later, the 13.5 s periodic pulsations were also detected in EXOSAT observations during 
non-flaring out-of-eclipse state (Pietsch et al. 1985). Nevertheless, while the periodic pulsations 
have been revealed several times during flaring events (Levine et al. 1991, Woo, Clark &; Levine 
1995, Woo et al. 1996), only once periodic pulsations have been observed during non-flaring state. 

The orbital period decreasing rate is -Porb/-Porb = (— 5.3=b2.7) x 10~^ yr~^. This value is smaller 
than that measured for Cen X-3 and SMC X-1 placing a lower limit on the rate of orbital energy 
dissipation by tidal forces. This fact, together with the relatively large radius of the companion 
(14 < R < 20 Rq), suggests that the companion is expanding rapidly, still burning hydrogen into 
a shell (Levine et al. 1993). In fact the resulting increase of the moment of inertia could acts 
to reduce the rotation rate of the companion star producing a relatively large difference between 
stellar rotation period and orbital period. This enhances the tidal friction effects producing a large 
orbital decay. 

Pulse-phase resolved spectroscopy from Ginga and ROSAT observations in the energy range 
0.2-30 keV shows that the phase dependent spectrum can be modeled as the sum of three com- 
ponent (Woo et al. 1996): a high energy component represented by a power law with photon 
index a ~ 0.7 with a high energy cutoff at 16.1 kcV and a folding energy of 35.6 kcV, a thermal 
bremsstrahlung component with temperature of 0.35 keV and a pulse-phase independent black 
body (kT ~ 0.03 keV). An iron line is also present at E ~ 6.6 keV. 

Several observations of the spin period evolution showed both spin up and spin down episodes. 
This suggests that the pulsar has reached a spin equilibrium state, wherein its spin period approx- 
imately equals the Keplerian period at the inner edge of the accretion disk. Because of the high 
X-ray luminosity of LMC X-4 and the relatively slow spin period, a very strong magnetic dipole 
has to be present (> 10^^ G cm^) (Woo et al. 1996; Naranan et al. 1985). 

Up to date no clear evidence of a cyclotron absorption feature has been reported. This is in 
line with the fact that the value of the cyclotron line energy deduced from the spin equilibrium 
(~ 100 keV) is out of the energy range (typically up to 30 kcV) of most X-ray satellites which 
antedate BeppoSAX. In fact the centroid of the line is related to the magnetic moment by the 
relation E'cyci = ^^-Q/J'So/ Rei^ + z) where usq is the magnetic dipole moment in units of lO^'' G 
cm^, Rq is the neutron star radius in units of 10^ cm, and (1 -|- 2;) = (1 — 2GMns/-Rc^)~^/^ is the 
gravitational redshift factor. Adopting a standard neutron star mass of 1.4 Mq with a standard 
radius of 10^ cm and fi > 10^^ G cm^, we obtain -Ecyci ^ 90 keV. The broad spectral band of the 
BeppoSAX satellite, extending up to 200 kcV, gives us the possibility of studying the spectrum up 
to the energies where the cyclotron absorption feature is expected. 
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In this paper we report the spectral analysis of LMC X-4 out of eclipse in high state observed 
by BeppoSAX in the 0.12-100 kcV energy range. We confirm the complex multicomponcnt nature 
of the spectrum of LMC X-4. Moreover we find that there are evidences of the presence of a 
cyclotron line, with the centroid around 100 keV, at a high degree of confidence. 

2. Observations and Light Curves 

BeppoSAX observed LMC X-4 with its narrow-field instruments (NFI; Boella et al. 1997a) 
in 1998 from October 20th 22:40 to October 22th 08:05 (UT) corresponding to the phase interval 
0.66-0.70 of the 30.3 days cycle, when the source was in its high state. The ephemeris (period: 
30.3 lb 0.05 days, epoch: 50087 MJD, corresponding to the beginning of the low-state phase) has 
been obtained by using the observations of LMC X-4 performed by the ASM on board RXTE. 

The BeppoSAX observatory covers more than three decades of energy, from 0.1 to 200 keV. 

The payload is composed by four coaligned instruments: the Low-Energy Concentrator Spectrom- 
eter (LEGS, 0.1-10 keV; Parmar et al. 1997), the Medium-Energy Concentrator Spectrometer 
consisting of two units (MECS, 1-10 keV; Boella et al. 1997b), the High-Pressure Gas Scintillation 
Proportional Counter (HPGSPC, 4-120 keV; Manzo et al. 1997), and the Phoswich Detection 
System (PDS, 15-300 keV; Prontera et al. 1997). The fields of view (FOV) of the instruments 
is approximately 1°. Since LEGS and MECS have imaging capabilities, we extracted the data 
from circular regions in the LEGS and MEGS FOVs of 8' and 4' radius, respectively, centered on 
the maximum of the point-spread function (PSF). The background subtraction is obtained using 
blank-sky observations and extracting the background spectra from a circular region corresponding 
to that one used for the source. The HPGSPC and PDS do not have imaging capabilities. In 
this case the subtraction of the background is obtained using off-source data collected during the 
rocking of the collimators. The eff'cctivc exposure during the observation arc so summarized: ~ 28 
ks for the LEGS, - 66 ks for the MEGS, ~ 29 ks for the HPGSPC, and ~ 39 ks for the PDS. 

In Figure 1 (upper panel) we show the light curve of LMC X-4, binned at 300 s, in the energy 
band 1.5-10.5 keV (MEGS data). The light curve shows that the X-ray source is eclipsed for 
about 20 ks. No important flaring episodes are observed in this light curve. The corresponding 
hardness ratio, i.e. the ratio of the counts in the 4.5-10.5 keV energy band to the counts in the 
1.5-4.5 keV energy band, is also shown in Figure 1 (lower panel). Although the source flux increases 
by ~ 20 — 30% during the pre-eclipse phase, the hardness ratio appears to be quite constant during 
our observation. Actually the hardness ratio decreases from ~ 2.5 to ~ 1 in a short time interval 
just after the eclipse. In the following analysis we considered only data out of eclipse, i.e. we 
excluded the time interval from 8.7^ s to 1.1 x 10^ s, containing both the eclipse and the short 
interval where the hardness ratio varies signiflcantly. 

The average flux out of eclipse in the high state is 9.7 x lO^^'^ ergs cm~^ s~^ in the 0.1-100 keV 
band. Adopting a distance to the source of 50 kpc, the resulting luminosity is ~ 3 x 10'^^ ergs s~^. 
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In this paper we concentrate our analysis on data out of eclipse. 

Because of the better statistics, we used MECS data to investigate for periodicities. We 
searched for periodic pulsations in a large interval (12.5-14.5 s) around the value of the spin period 

reported in the literature, 13.5 s, by using folding techniques. No pulsations have been found. In 
fact the maximum excess that we have found had = 93.82 for 63 degrees of freedom, which is 
not significant. 



Spectral analysis was performed on the average photon spectrum of LMC X-4, out of eclipse, 
in the energy range 0.12-100 keV. The energy ranges used for each NFI are: 0.12-3 keV for LEGS, 
1.8-10 keV for MECS, 9-30 keV for HPGSPC, and 15-100 keV for PDS. Different normalizing 
factors for the four instruments were included, fixed to 1 for the MECS and kept free for the other 
instruments. 

The continuum is described by a power law with a high energy cutoff corrected by photo- 
electric absorption. Indeed, by using the XSPEC models, POWERLAW and HIGHECUT, an 
evident residual remains at the cutoff energy. Since it just appears at the cutoff energy, where 
the fitting function shows a cusp, this residual probably is an artifact of the model. To verify this 
hypothesis we changed the XSPEC HIGHECUT function, whose mathematical form is: 



with a function F{E) that avoids the cusp. In fact F{E) has a smoothed region of width W, 
between the low energy range (E < EcutoS — W/2) and the high energy range {E > -BcutoflF + W/2) 
described by a third degree polynomial function, in order to obtain a function that is continuous 
with its derivatives (sec also Burderi et al. 2000); The mathematical form of our function is: 



where A, B, C, D are constants calculated imposing continuity conditions for the function F(E) 
and its derivatives while W is left as a free parameter of the fit. Using this model, the residuals at 
the cutoff energy disappear, and the fit improves significantly. 

An evident soft excess remains at low energies below ~ 2 keV. A simple black body model 
is not a good description of the soft excess, giving x^/d.o./. ~ 727/537. Following the literature 
(Woo et al. 1996) we fit this feature with a thermal bremsstrahlung emission plus a black body. 
The resulting parameters describing the model are reported in Table 1 (model 3 and 4) and can be 
so summarized. The photon index of the power law is ~ 0.6, the cutoff energy is ~ 18 keV and the 



3. Spectral Analysis 




for E < EcutoS 
for E > -Bcutoff 



(1) 



' 1 for ^ < £;eutoff - VF/2 

F{E) = 1 AE^ + BE^ + CE + D for S,,tofr -W/2<E< Eeutoff + W/2 
exp{ %2l-f ) for E > ^cutoff + W/2 



(2) 
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folding energy is ~ 17.5 keV. The photo-electric absorption corresponds to an hydrogen column of 
~ 5 X 10^*^ cm~^. The temperature associated to the bremsstrahlung is 0.85 keV while that of the 
black body is ~ 6 x 10~^ keV. 

An alternative way of fitting the soft excess is to model it with a Comptonization model (see 
Table 1, model 1 and 2) . We used the COMPST model of XSPEC (v.lO). We preferred COMPST 
to COMPTT because COMPST needs less free parameters. In addition, using COMPTT, the 
temperature of the seed photons obtained by the fit is less than 0.1 keV, i.e. it is out of the 
spectral range, making this estimate unreliable. The temperature associated with the hot electrons 
is ~ 1 keV and the relative optical depth is ~ 16. The parameters of the other model components 
remain substantially unchanged. 

The presence of a gaussian emission line at ~ 6.5 keV is evident. Indeed the reduces 
significantly if we model it with two gaussians centered at the energies 6.1 and 6.6 keV (equivalent 
width ~ 200 and 65 cV, respectively). Comparing the model with one gaussian and that with two 
gaussians, the F-test gives a probability greater than 0.997 that the improvement is not casual. 
The simultaneous presence of a narrow line at 6.6 keV and a broader one at 6.1 keV implies two 
different iron ionization stages. In addition an absorption edge is observed at about 1.25 keV with a 
very large statistical evidence (the F-test probability of chance improvement is less than 1 x 10~^). 
This feature could be attributed to the absorption edge of Fe XVII (1266 eV). 

Large residuals also appear in the hard part of the LMC X-4 spectrum starting from 40 keV 
and extending up to 100 keV (see figure 2, upper panel). We interpreted this feature as a cyclotron 
absorption line and modeled it with a gaussian profile, i.e. we used a multiplicative function given 
by: f{E) = 1 — ^cyci exp[—{E — £^cyci)^/(2c"^)], with the condition that f{E) = when it assumes 
negative values. Since the width of cyclotron lines is most probably of thermal origin (see e.g. 
Mihara 1995; Dal Fiume et al. 1999; Burderi et al. 2000), we preferred to use a gaussian profile, 
instead of a lorentzian profile, to fit this feature. Using this model the residuals disappear (see 
Figure 2, lower panel). The centroid energy of the line is ~ 100 keV and the sigma is ^ 50 keV. 
The addition of this component reduces the significantly, as it is evident comparing models 1 
and 3 with respect to models 2 and 4 in Table 1. The F-test gives a probability of 4 x 10~^ for a 
chance improvement of the fit. 

Summarizing the whole 0.1-100 keV spectrum of LMC X-4 can be described as follows: 

S{E) = WABS X EDGE x CYCL x {SOFT + 2 GAUSS + F{E) x POW) (3) 

where WABS is the photo-electric absorption, EDGE is the absorption edge, CYCL is the cyclotron 
absorption feature modeled as a gaussian, SOFT is the component that can be modeled as a 
bremsstrahlung plus a black body or as a Comptonization, 2 GAUSS represents the two iron 

emission features at about 6.5 keV, F(E) is the cutoff function above described and finally POW 
is the power law component. The fit corresponding to the model with bremsstrahlung and relative 
residuals are plotted in Figure 3, while the fit corresponding to the model with Comptonization 
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and relative residuals are plotted in Figure 5. In both cases the cyclotron absorption line is present. 
The unfolded spectra corresponding to the two scenarios are drawn in Figure 4 and 6, respectively. 

4. Discussion 

We have found evidence, with a high degree of confidence, of a high energy cyclotron absorption 
feature in the broad band (0.1-100 keV) spectrum of LMC X-4. Modeling this feature as a gaussian 
in absorption, we have found the centroid energy at ~ 100 keV with sigma ~ 50 keV, corresponding 
to a magnetic moment of 1.1 x 10^^ Gauss cm^. The addition of this new component to our model 
significantly reduced the x^, with a probability of chance improvement of 4 x 10~^. This is in line 
with the expectation that LMC X-4 should have a very strong magnetic field. 

There are some evidences that LMC X-4 have reached its spin equilibrium wherein its spin 
period is close to the Keplerian period at the inner edge of the accretion disk, which is probably the 
magnetospheric radius of the system. In fact several observations show both spin-up and spin-down 
regimes (see e.g. Woo et al. 1996). The average values of -Ppuise range from —4. Ox 10~^ s/yr between 
the HEAO 2 and the first of the EXOSAT observations to +1.9 x 10"^ s/yr between EXOSAT and 
the most recent Ginga and Rosat observations. The radius i?eq at which the Keplerian period of the 
orbiting matter in the disk equals the pulsation period Ppuise is simply i?eq ~ 1.5 x ^^^fn^^^P^Use 
cm, where m is the neutron star mass in units of solar masses and Ppuise is the spin period in 
seconds. For Ppuisc = 13.5 s, Pcq results ~ 9 x 10^ cm. In addition we know that, for an accreting 
magnetized neutron star, the magnetosphere radius, Pm, is (see, e.g., Burderi et al. 1998): 

Pm ~ 4.3 X L'^^'' m^/^ Rq'^''' cm (4) 

where is the magnetic moment in units of lO^'^ Gauss cm^, L37 is the intrinsic luminosity in units 
of 10^^ erg/s, Rq is the neutron star radius in units of 10^ cm, and finally (p (< 1) is the correction 
to the Alfven radius for the presence of the accretion disk. Equaling Peq and Pm it is possible to 
estimate the strength of the magnetic moment. Taking into account the intrinsic high luminosity of 
LMC X-4 (~ 3 X 10^^ crgs/s in the 0.1-100 keV band), we derive a high magnetic dipole moment, 
of the order of 2 x 10'^^ Gauss cm^. If this is correct, any cyclotron spectral feature in this source 
has to be expected at energies around 150 keV. Indeed we find evidences of an absorption feature 
in the hardest part of spectrum of LMC X-4 centered at 100 keV which is consistent with the 
expected value. 

The width (gaussian a) of the cyclotron line that we found is ~ 57 keV. If the broadening is 
caused by thermal Doppler effects, the relationship among the folding energy of the exponential 
cutoff, -Bfoid) that is a measure of the temperature close to the neutron star, the centroid energy of 
the cyclotron line, -Ecych and its width, a, is: a ~ Pcyci(Pfoid/''TT'eC^)^''^- Adopting this relationship 
the expected width of the line would be ~ 26 kcV that is almost half the measured broadening of 
the line. However, note that the width of the cyclotron line might be not well constrained since 
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the centroid energy of the line is at the end of the energy range of PDS, and therefore only the left 
wing of the line is visible. 

The whole 0.12-100 keV spectrum of LMC X-4 is complex. The primary continuum is well 

fitted by a power law with a high energy cutoff, smoothed at the cutoff energy. This component 
could be produced by Comptonization of the radiation coming from the hot spots on the neutron 
star by the plasma in the accretion column, very close to the neutron star surface. Anyway, at 
low energies, an evident soft excess remains. We fitted this excess with two different models: 1) 
bremsstrahlung emission with the addition of a black body, 2) Comptonization. 

In the first scenario, the black body component is probably emitted by the accretion disk, 
whereas the bremsstrahlung component comes from an optically thin plasma surrounding the sys- 
tem. In fact, considering that the luminosity and the temperature of the black body component, 
obtained by the fit, are 1.3 x 10^^ ergs s~^ and 6.3 x 10~^ kcV, respectively, the typical radius 
associated with emitting region of the black body component is ~ 2.5 x 10^ cm. This value is in 
agreement with that found for the magnetospheric radius under the spin equilibrium hypothesis. 
With regard to the bremsstrahlung component, the calculated luminosity is 2.2 x 10^^ ergs s~^. 
This give us the possibility to estimate the radius of the bremsstrahlung emitting region. The 
luminosity of the bremsstrahlung component, i^brems) can be expressed as: 

i^brcms ~ 5.8 X 10-24ri/2iVf y erg/s (5) 

where T is the electron temperature in keV, A^e is the electron density and V is the volume of 
the emission region in cm^. Considering the relationship between the optical depth r and iVg: 
T ~ (7T-/Ve-Rj where (Jt is the Thomson cross section and R is the radius of the emission region, and 
supposing a spherical geometry, we obtain: 

ibrems ~ 5.5 X IO^^T^/^Rt^ crg/s (6) 

Adopting the electron temperature of 0.86 keV obtained by the fit, the radius of the emission region 
IS 4.9 X lO^V^^ 

cm. For an optical depth associated to the bremsstrahlung not greater than 1, 
the resulting radius of the emitting region is larger than the radius of the Roche lobe containing the 
neutron star. Therefore, in this scenario the black body component could be emitted by optically 
thick plasma at the magnetosphere, and the bremsstrahlung component could be produced by 
plasma of the stellar wind from the companion star, ionized by the X-ray emission from the pulsar. 

It is also possible to model the soft excess with a Comptonization model. The corresponding 
luminosity of this component is 5 x 10^^ ergs/s. Knowing the temperature of the seed photons 
and that of the hot electrons, the luminosity of the emitting region and its optical depth, it would 
be possible to infer the characteristic radius of the seed photons emitting region, where the seed 
photons originate. The COMPST model does not give the possibility to obtain the temperature 
of the seed photons from the fit, therefore we derived this temperature assuming that it is the 
temperature at the inner radius of the accretion disk (assuming the standard model of Shakura 
Sz Sunyaev 1973), under the hypothesis that the disk is truncated at the magnetospheric radius. 
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Using the value of the magnetic field corresponding to the cyclotron line, we obtain a temperature 
of the seed photons of 1.8 x 10~^ keV. From this we obtain that the radius of the region emitting 
the seed photons is 36 x 10^ cm. This value is in agreement within a factor 5 with that obtained 
from equation 4 where we have adopted = 11, as deduced by the fit of the cyclotron feature. 
However, note this model might be unadequate to describe the whole emission mechanism, since 
at this cold temperatures and large optical depths, free-free absorption, that is not included in the 
COMPST model, could give a not negligible contribution. 

An iron emission line is detected at about 6.5 keV. we find that it can be fitted by two 
gaussians at energies 6.1 and 6.6 keV, respectively. The gaussian at lower energy has a large sigma 
(~ 0.95 keV), whereas the gaussian at higher energies is much narrower. We note that the width of 
the broad line component cannot be explained by the Doppler effect induced by the fast rotation 
at the magnetospheric radius. In fact considering a magnetic moment of 1.1 x 10'^^ Gauss cm^ 
(as estimated by the cyclotron line energy) we obtain Ai/ ~ 0.1 keV at 6.4 keV, that is much 
smaller than the measured width. Complex iron line features are not unusual in high mass X-ray 
binaries (see e.g. Audley et al. 1996, Ebisawa et al. 1996). The energy values that we obtained 
are consistent with the iron features expected at 6.4 and at 6.7 keV, within the associated errors. 
The usual interpretation is that the line at 6.7 keV is produced by highly ionized iron and probably 
is emitted in a hot corona, while the line at 6.4 keV corresponds to neutral iron and might come 
from the radiation reflected by the accretion disk. In this case we would expect the 6.7 keV line to 
be broad because of the Compton smearing in the hot corona, and the 6.4 keV line to be narrow 
given that the disk is truncated at a large radius by the presence of the magnctosphere. We observe 
that our data of LMC X-4 indicate an opposite behaviour, with a broad 6.4 keV line and a narrow 
6.7 keV line. Probably this feature can be explained by line blending, for the presence of several 
iron line components with energies between 6.4 and 7 keV, or by the presence of a line with a 
narrow core but broad wings. Certainly other high resolution observations in the iron line region 
are needed to address this question. 

All the models reported above need an absorption edge at about 1.25 keV. The nature of this 
feature is highly uncertain. The energy of the edge is close enough to the K-edge of high ionized 
Ne or neutral Mg, or to an L-cdgc of moderately ionized Fc. If we attribute this feature to FeXVII 
(edge at 1.266 keV), the optical depth corresponding to this feature is: 



where <;;iwith-edge ^■iid ^no-edge ^re the flux with and without the absorption edge, respectively, 
calculated in the 1.25-3.0 keV range. Assuming that the cross section of this process is comparable 
with that of the ejection of electrons from the K-shcU [a --^ 1.7 x 10~^^ cm^) the corresponding iron 
column is ~ 0.6 x 10^^ cm~^. This value could be compared with that deduced from the hydrogen 
column obtained by the fit of the photoelectric absorption. In that case the hydrogen column is 
~ 5 x 10^° cm~^ so the deduced iron column (assuming solar abundances) could be of the order of 
5 X 10^^ cm~^. The good agreement between the two values indicates that the edge at 1.252 keV 




(7) 
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could be a real feature produced by the iron around the neutron star. Other authors have also 

reported similar features. Parmar et al. (2000) have found an absorption edge in the spectrum of 
the low mass X-ray binary X 1822-371 at about 1.33 keV and optical depth 0.26. Moreover laria 
et al. (2000) have observed the edge in the spectrum of the source 4U 1254-69 at 1.26 keV and 
optical depth 0.15. 

This work was supported by the Italian Space Agency (ASI), by the Ministero della Ricerca 
Scientifica e Tecnologica (MURST). 
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Table 1: Fit of LMC X-4 spectrum during the high state in the energy band 0.12-100 keV. Un- 
certainties are at 90% confidence level for a single parameter. Tedge is the absorption depth of 
the edge at the threshold. kTee is the temperature of the black body. The normalization of the 
black body, Nbb is in units of Lsg/DgQ where L36 is the luminosity of the source in units of 1036 
crgs/s and D50 is the distance of the source in units of 50 kpc. kTBrcms is the temperature of the 
plasma emitting for bremsstrahlung. The normalization of the bremsstrahlung, NBrems ) is in units 
of (1.27 X 10-^^/(47ri:'|o)) / ^e^idV, where Ne and Nr are the electron and ion densities (cm'^), 
respectively. kTcompst is the temperature of the hot electrons Comptonizing cool photons. The 
COMPST normalization is in XSPEC (v. 10) units. The power law normalization is in units of 
photons keV"^ cm"^ s"^ at 1 keV. The (emission lines) gaussian normalization is in units photons 

— 2 —1 

cm s . 



Parameters 


model 1 


model 2 


model 3 


model 4 




O.U_o,4 


5.0 ±0.5 


5± 1 


5±1 


Edge Energy (keV) 

'^edge 


1.28 ±0.04 
0.4 ±0.1 


1 90+0.02 

-'■•^°-0.03 

0.4 ±0.1 


1.25 ±0.04 


1.25 ±0.04 
0.5 ±0.1 


K i BB {KG V ) 

Nbb 

kTBrcms (keV) 

NBrems 






[p.o It i.O j X iU 
2.4+t^ X 10-2 

0.81^:3 

Atl X 10-2 


1^0. =t i.OJ X iU 

2.4t4-2 X 10-2 
0.8 ±0.2 

(4± 1) X 10-2 


kTeoxnpst (keV) 

'^compst 
■f^compst 


q+4 

16 ±5 
2.ltoi X 10-2 


9+1° 
c+e 

1.7 ±0.2 


- 


- 


Photon Index 

Npower law 

Smooth Width (keV) 
Scut (kcV) 
Efoid (keV) 


U.O_g 2 

(1.1 ±0.2) X 10-2 
9±4 

18 ± 1 
18 ±3 


0.5 ±0.1 
7;2 X 10-3 

17 ±3 

18 ± 1 
12.2 ±0.4 


0.6 ±0.1 

i.ol°:2 X 10-2 

9±4 
18 ±1 


64+" °3 

"■"^-0.04 

9.9to l X 10^ 
15 ±3 
18.7±0.5 
12.7 ±0.3 


Epe (keV) 1st gaussian 

CTpo (keV) 

Jpe (ph cm-2 g-i) 

EWpe (eV) 

Epe (keV) 2nd gaussian 

crpc (kcV) 

Ipe (ph cm-2 s~^) 
EWpe (eV) 


6.1 ±0.3 

u.y_Q 2 

9tr X 10-' 

169 
6.6 ±0.1 

(0.1 ±0.2) X 10-2 
3^f X 10-4 
64 


1 1+0.5 

-■^•^^-0.2 

8tl X 10-'' 
183 
6.6 ±0.1 

(0.1 ±0.1) X 10-2 
(2 ± 1) X 10-* 
64 


6 2+°-2 

"•^-0.8 

0.8 ±0.4 
6tf X 10-4 
215 
6.6 ±0.1 
7tf X 10-2 
2.3lo.4 X 10-4 
67 


6 2+°-2 
"•^-0.7 

6tl X 10-4 
267 
6.6 ±0.1 
7+f X 10-2 
2+1 X 10-4 
69 


i^cyci (keV) 
CTcyci (keV) 

Acycl 


lootr 

60ir 

i.ot^l 




4^+140 

1 n+°-4 

'-■^-0.3 






614/536 


636/539 


614/535 


638/538 
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FIGURE CAPTIONS 

Figure 1: Upper panel: LMC X-4 light curve in the energy band 1.5-10.5 keV band (MECS data). 
The bin time is 300 s. Lower panel: Corresponding hardness ratio [4.5-10.5 keV]/[1.5-4.5 kcV]. 
Figure 2: Upper Panel: Residuals in unit of a with respect to model 2 of Table 1 (without the 
cyclotron absorption line correction). Lower Panel: Residuals in unit of cr with respect to model 
1 of Table 1. In this case the cyclotron absorption line correction is taken into account. Different 
symbols have been used to distinguish between HPGSPC data (filled triangles) and PDS data (open 
triangles). 

Figure 3: Pulsc-avcraged count spectrum (0.12-100 kcV) of LMC X-4 and residuals in units of 
cr corresponding to model 3 of Table 1. In this case the soft part of the spectrum is fitted with 
brcmsstrahlung plus black body. The cyclotron absorption line is also included. 
Figure 4: Unfolded energy spectrum and best fit model corresponding to model 3 of Table 1. The 
solid line with the data on top is the total spectrum, the dashed lines are the blackbody (left) and 
the bremsstrahlung component (right) respectively, the dot-dot-dot-dashed line is the power-law 
component, the dot-dashed line is the gaussian at 6.1 keV, and the dotted line is the gaussian at 
6.6 keV. 

Figure 5: Pulse-averaged count spectrum (0.12-100 keV) of LMC X-4 and residuals in units of 
a corresponding to model 1 of Table 1. In this case the soft part of the spectrum is fitted with a 
Comptonization. The cyclotron absorption line is included. 

Figure 6: Unfolded energy spectrum and best fit model corresponding to model 1 of Table 1. The 
solid line with the data on top is the total spectrum, the dashed line is the comptonized emission 
component, the dot-dot-dot-dashcd line is the power-law component, the dot-dashed line is the 
gaussian at 6.1 keV, and the dotted line is the gaussian at 6.6 keV. 
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